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1 Introduction 

The idea of a chiral spin liquid state spontaneously breaking parity (P) and time reversal (T) invari- 
ance has attracted considerable interest recently. It has first been proposed as a possible ground state 
for the two dimensional •S' = ^ Heisenberg model on a square lattice frustrated with a sufficiently 
strong antiferromagnetic next-nearest neighbour interaction [|1]]. Subsequent studies of this model 
have found an enhancement of a chiral order parameter, comparison with other possible states how- 
ever suggests that the chiral spin state is unstable [^. Different lattices, in particular the triangular 
and Kagome one, have also been studied, however no firm evidence of a chiral spin state has been 
found so far. Here the frustration is a consequence of the lattice geometry. 

To characterize a chiral phase several 'order parameters' have been introduced: For lattices built 
from triangular plaquettes the vector chirality [^ 

X = SiX S2 + S2XS3 + S3XS1 (1.1) 

with Sj being the three spins on the comers of a triangular cell has been discussed. While a sponta- 
neous symmetry breaking in X appears to be unlikely in an isotropic Heisenberg system its properties 
have been studied in stacked triangular antiferromagnets with XXZ type anisotropy [^ ^ which are 
realized in the ABXs-type compounds such as CsCuCla (see articles in [||]). 

A rotationally invariant operator which has a non zero expectation value in a phase with broken 
P and T symmetry is defined through [|1]] 

X(,,3> = Si ■ (S2 X S3) . (1.2) 

Finally, a topological "Chems" number measuring the dependence of a quantum state on a twist in the 
boundary conditions has recently been introduced by Haldane and Arovas and used to characterize 
the ground state of a Heisenberg model on a hexagonal lattice subject to PT breaking fields [0]. The 
actual computation of the Cherns number, however, is restricted to rather small systems thus limiting 
its use in studies of a phase diagram at present. 

While the existence of a phase with spontaneous broken chirality in a frustrated two-dimensional 
Heisenberg model has not been established so far (a possible candidate may be the Kagome lattice 
[^) the experimental studies of ABX3 compounds indicate that the frustration leads to a rich phase 
diagram if a magnetic field is applied [||, ^ |TOp . 

Lacking a model with spontaneous broken PT-symmetry it is useful to consider models con- 
taining terms breaking these symmetries explicitly. Studies of such systems allow to gain a better 
understanding of the properties of the chiral spin liquid state and means for its characterization. In 



this paper we analyze a system of two spin-| Heisenberg chains coupled by exchange terms on diag- 
onal bonds as shown in Fig. |1] described by the Hamiltonian 

2N 

'Ho = 2_^ 1^1 Sn ■ S„+i + J2^n ■ S„+2} (1-3) 

n=l 

with periodic boundary conditions S2Ar+fc = S^ and even and odd indices labeling spins on the two 
sub chains respectively. For J2 = this is just the Bethe Ansatz soluble Heisenberg chain [fTT|]. As 
long J2 < J2C ~ 0.25 Ji the model continues to have gapless excitations above a translationally invari- 
ant ground state just as the single chain. Increasing the frustrating next-nearest neighbour interaction 
beyond J2C the system has two degenerate dimerized ground states leading to the Majumdar-Ghosh 
model with a simple dimer configuration for its ground states at J2 = \Ji [0. 



To force the system into a chiral spin state we add terms breaking PT-symmetry explicitly 



1 ^ 

^v = - ^ {Xl S2n-1 ■ (S2n X S2n+l) " X2 82™ " (S2n+1 X S2n+2)} ■ (1-4) 

n=l 



'•^ 2 



Such multispin exchange terms may in fact be relevant for the description of certain realizations of 
the two dimensional Heisenberg model such as ^He layers on a graphite substrate [[T3p. Note, that 



the invariance under translations by one lattice site {n —^ n + 1) h destroyed by these term unless 
Xi = — X2- For the full Hamiltonian including the chiral terms one has additional integrable models: 
for a "staggered" chiral field xi = — X2 the operator H^ is one of the hierarchy of integrals of motion 
of the Heisenberg chain, thus commutes with 7io(-^2 = 0). The consequences of the competition 
between these operators have been studied in [jT^. Choosing the parameters in 

n = no + n^ (i.5) 

as 



Ji = 2(l-fi:), J2 = ^, Xi = X2 = 2VM1 - /€) (1.6) 

one obtains a family of integrable models of generalized spin ladders introduced recently [|15], [1^, |17p. 
Varying of the free parameter k from to 1 the system evolves from a single Heisenberg chain to a 
pair of decoupled ones, changing the sign of k reverses the sign of the chiral field, but doesn't affect 
most other properties of the system. 

Here we investigate the properties of the ground state and excitations of this model. In the follow- 
ing section the ground state energy and spectrum of the low lying excitations of the integrable model 
( |1.6[ ) are computed exactly. As in the Heisenberg chain the excitations over the antiferromagnetic 
(singlet) ground state are found to be spinous coming in pairs. For the characterization of the chiral 



properties we note that due to its quasi one dimensional character it is not possible to define an ana- 
logue of the topological Cherns number for this system. In the following we choose the expectation 
value of a uniform extension of the chirality ( [L2| ) as a measure of the chirality 

1 ^ 

X = T7 2_^ {S2n-1 • (S2n X S2n+l) — S2n, " (S2n+1 X 8271+2)} • (1-7) 

n=l 

Unfortunately, expectation values of operators not commuting with the Hamiltonian such as x are not 
easily accessible within the framework of the Bethe Ansatz. Our results regarding (£) are obtained 
from numerical diagonalization of finite clusters. In Section ^ we study the phase diagram of the 
integrable model subject to a uniform external magnetic field h || 5. A characterization of the phases 
based on counting the number of gapless excitations supported by the system is possible from the 
Bethe Ansatz analysis. We identify three different phases in the n-h plane: for sufficiently large 
h > hc2 the system shows saturated ferromagnetism. For ^ < k < 1 a phase with low lying 
excitations at four different wave numbers ±A;i 2 is found for magnetic fields hd < h < hc2 (a 



similar phase diagram has recently been established QISQ in an integrable chain of alternating spins 
I and 1 [|T9p). To characterize these phases we can compute the magnetization again from the Bethe 
Ansatz while we have to rely on numerical results from finite systems for the chirality. In the final 
section we summarize our findings and comment on the properties of systems where the parameters 
Xi 2 are tuned away from the integrable point. 

2 Ground state and excitations of the integrable model 

As mentioned above, choosing the exchange constants for the Hamiltonian ( |1.5p as in ( |1.6D gives a 
system which is integrable by Bethe Ansatz methods. Starting from the ferromagnetic state with all 
2N spins pointing up one can reduce the solution of the Schrodinger equation in the sector with M 
overturned spins to a system of algebraic equations 

for the M complex rapidities \j. Each solution of these Bethe Ansatz equations ( |2.1D corresponds 
to an eigenstate of H with spin S = N — M = S^. Up to an overall constant the corresponding 
eigenvalues are given by 

E{M) = E (^^o(A, + |) + eo(A, - |) j + Mh, e~o(A,) = ~^^^^ ■ (2-2) 



Here we have used the fact, that we are constructing eigenstates for fixed S^, to include the effect of 
an external magnetic field in the Hamiltonian H ^ H — hS^. 



A generic solution {Xj} of ( |2.1D is organized groups of uniformly spaced complex rapidities, so 

called strings 

,(m) ,. m + lm + 3 m-1 

X)'=x + tfij, Hj = —, ^'•••'^~- *^2-3^ 

In the thermodynamic limit the ground-state is made up of real A's only (1-strings). Their density 
p(A) is given in terms of a linear integral equation 

^ 1^1 1 



P(A) + J_^ dM-(A -,)M-^ I (ITlpTI ^ (I^lpTT j ■ <^-^' 

The kernel of this integral equation is K{\) = ^ pxi- The dependence on the magnetic field is 
incorporated in the value of the integration boundaries A which are to be chosen such that the total 
density is J_^ d\p{\) = (M/N). For a vanishing external magnetic field h = the ground state can 
be shown to be a singlet (M = N) and the rapidities fill the entire real axis resulting in A = oo in 



([2.4D. Hence p can be computed by Fourier transform, resulting in 

,(A) = p(AH-|)+p(A-§), „ithp(A) = ^^^. (2.5) 

From ( p^ the ground state energy Eq per spin is given by [ [T7| ] 

^ = -y^rfAp(A)6o(A) = y^rfA(p(A + -)eo(A) + (p(A--)6o(A)) 

= -\n2-\mi-'^K)-^{^--'-~K) + ni+'^f^)-nl+'f))- (2.6) 

The term containing digamma-functions \l'(x) increases from — In 2 to as k varies between and 
1. 

Low lying excitations of the system are parametrized by holes in the distribution of real {Xj}- The 
dispersion of these spinous is determined by the integral equation for the dressed energies: 

In this grand canonical approach the ground state is characterized as the one in with all states with 
negative dressed energy e(A) are filled. For vanishing or sufficiently small (see below) magnetic fields 
this condition is related to that used in ( |2!4| ) above through e(±A) = 0. Again, we can solve ( p7| ) by 
Fourier transform for /i = 0, giving 

e(A) = e(A + J) + e(A - J) , e(A) "" 



2 cosh(A) 



A;(A) = ~k{\ + -) + ~k{\ - -) , k{X) = arctan(A) - - . (2.8) 



Eliminating A from these equations one obtains the spinon dispersion e{k) which is shown for several 
values of k, in Figure 0. 

Comparing these results with the known solution of the usual S = ^ XXX Heisenberg chain [ ^ 



we find that they coincide in the limits k -^ 0, 1 as is expected from the Hamiltonian. Zero tem- 
perature quantities such as densities of rapidities or their dressed energies are simply superpositions 
of the corresponding quantities for the XXX chain with argument shifted by k (see (|2.4D, (|2.7|)). On 
the basis of the properties of the low lying excitations this comparison can be extended to the critical 
properties of the system: For zero h the continuum limit of the system can be identified as a con- 
formal field theory (CFT) with central charge c = 1 for any < k < 1. As for the XXX chain we 
expect this situation to be described by a level- 1 SU(2) Wess-Zumino-Witten model. For k = 1 there 
appears another massless mode leading to low energy properties corresponding to two c = 1 models 
(see also Section^ below). 

Finally we have studied the effect of the symmetry breaking terms 7i^ in the integrable Hamil- 
tonian. It is instructive to start with the exact solution of a system of four spins (see also [|l]]). The 
Hamiltonian is given by (lO) with N = 2 and the couplings are 



Ji = l, J2, Xi = X2=:«:, <«:<!. (2.9) 

From the exact solution we know that the ground state of the integrable model is always a singlet. 
The same holds for the model given by (|Z9|). Hence it is sufficient to diagonalize the Hamiltonian in 
the two-dimensional singlet subspace. Depending on J2 there are two distinct cases: 

1. J2 = ^ ■ At the Majumdar-Ghosh point the singlets are degenerate for vanishing chiral field 
whereas for finite k the operator £ ( p77| ) lifts this degeneracy. One also observes that the Hamiltonian 
simplifies to H = S^ + kx (up to a constant) and [S^,x] = 0. Therefore the chirality x can be 
diagonalized in the singlet subspace, yielding a ^-independent chirality which is found to be one of 
its eigenvalues (x) = ±|v^. 

2. J2 7^ I : Tuning J2 away from ^ the ground state of the k = 0-system is unique. At this point the 
Hamiltonian is P- and T-symmetric leading to a vanishing chirality in the ground state. Switching 
on the chiral field the two singlet eigenstates of the Hamiltonian have non-zero expectation values 

(X): 

(2) = ±g ^^ ^ (2.10) 

^ ' 2 ^3fi:2 + 2(2J2 - ly 

The "chiral susceptibility" d{x) /Ok diverges for k = at the MG point. For k ^ 1 the chirality 
approaches its eigenvalue ^V^. 

For larger systems no exact results on the chirality can be obtained (as mentioned in the introduc- 
tion it is not possible to compute this expectation value directly from the Bethe Ansatz solution). For 



small systems (up to 24 spins) we have used a Lanzcos algorithm to compute some low lying states 
numerically. In Fig. ^ we present our results on ground state chirality from these data: (%) vanishes 
for K —> 0, 1 as is expected from the single chain Hamiltonian in these limits. For intermediate values 
of K we find that the finite size corrections to the chirality are very small, so that the properties of 
the infinite system can easily be read off from the numerical data. We find that the maximum of the 
chirality is obtained at k ^ 0.58 with about two thirds of the largest possible value v^/2. 

3 Magnetic phase diagram of the integrable model 



While the zero field properties of the model (|1.5|) are closely related to those of the single XXX 



Heisenberg chain, addition of an magnetic field gives rise to an interesting phase diagram (see Fig. 0): 
For small values of k we find that the magnetic field first breaks the SU(2) symmetry giving a c = 1 
Gaussian CFT with anomalous dimensions depending on the magnetic field. Increasing the magnetic 
field h beyond a value hdin) two additional "Fermi points" with gapless excitations arise, leading to 
a low energy spectrum as described by two c = 1 Gaussian CFTs. Finally for magnetic field h > hc2 
the ground state of the system saturates ferromagnetically. 

To calculate the corresponding phase boundaries as a function of k and the expectation values of 
the magnetization a = {S^) /{2N) as a function of the magnetic field we have to solve the integral 
equation (|2.7[). The numerical results for the magnetization are given in Figure]^. 

Analytical results can be obtained near the critical field hc2 which is determined from the condition 
that e(A) > for all values of A. For fields h ^ hc2 the integral equation can be solved by iteration 
which allows to study the nature of low lying excitations and the dependence of the magnetization on 
h. Here one has to distinguish three different cases: 

(i) For < K < J the minimum of the bare dispersion given by the driving terms in ( |2.7D is at 
A = (see Fig. ||(a)). The dressed energies are nonnegative for magnetic fields 

h>hc2 = 4(1 - k) . (3.1) 

In this region the magnetization reaches its maximum according to a square root law 



(T = l--J-r^Vhc2-h. (3.2) 

The system has massless excitations near the pseudo Fermi points e(±A) = 0. Letting h ^ one 
finds that there is no additional phase transition, hence we can identify these excitations with the 
spinons forming the excitation spectrum at zero magnetic field. Standard techniques [ Pj ] can be used 



to compute the spectra of finite size systems from the Bethe Ansatz equations (P?T|). The low lying 
energies compared to the ground state energy of the infinite system are (L = 2N) 

AE = -^v+'^v(A+ + A-) (3.3) 

bL L 

where v is the velocity of the massless magnons and the conformal dimensions of primary operators 

are 



ijlAMieADJ . (3.4) 



2 
AM is an integer denoting the change in 5^ induced by the operator, AD is an integer or half integer 

proportional to the momentum of the excited state (due to backscattering). The dressed charge ^ = 

^(A) is given in terms of the linear integral equation ^(A) + / ^ dfiK{\ — n)^{n) = 1. Depending 

on the external magnetic field it varies between 1/a/2 for /i = and 1 for /i — » hc2- The spectrum 



( |3.3D allows to identify the central charge c = 1 and the operator content of a Gaussian CFT with 
U(l) symmetry as mentioned above. 

(ii) K = J. proceeding as in (i) one finds hc2 = 'i- Unlike for smaller k the dependence of the bare 
dispersion on the spectral parameter near the minimum is oc A"^. This leads to a different dependence 
of the magnetization on the field near hc2'- 

a = l--{'i{K2-h))-^ . (3.5) 

TT 

The classification of the low lying excitations as well as the interpretation concerning the conformal 
properties of the system coincide with those of case (i) above. 

(iii) In the region | < /t < 1 two degenerate minima of the bare dispersion exist at nonzero 
values A = ±A*^°) of the spectral parameter (Fig. ^b)). This does not affect the dependence of the 
magnetization on the magnetic field for fields near the critical value 

hc2 = l + K-^ (3.6) 

As for K < i we have a square root dependence of a on /i ^ hc2 

a = 1 — const. \/hc2 — h (3.7) 

which a /t-dependent constant. The spectrum of low energy excitations, however, is different in this 
regime. For h < hc2 there are two filled "Fermi seas" [— A2, — Ai] and [Ai, A2] of quasiparticles near 
±A'^°) giving rise to two branches of massless spin excitations with different magnon velocities in the 
system. This situation is very similar to the one observed in the system with XY-type anisotropy and 



large staggered chiral field [inp: the expressions for the leading finite size corrections to the energies 



are (for the generalization of pl| ] to the case of several branches of excitations see also [^ p3| |) 

AE = -^ {v^ +V2) + ^- (t;i(A+ + Ar) + t;2(A+ + A^-)) . (3.8) 

bL L 



Here Vi are the velocities of the excitations (linearized near the Fermi points ±Aj) and the primary 
conformal dimensions A^ are found to be 

The numbers AM^ are the difference between the number of quasi particles in the Fermi seas in 
the ground state and excited state, respectively. AD^ counts the backscattering events in the excited 
state (giving rise to excitations with momenta being multiples of twice the Fermi momenta of the 
quasiparticles). All conformal dimensions can be parametrized by the four numbers ^f" = C^(Ai) 
which are given in terms of the integral equations 

e+(A) = 1 - / ' K{X - /.)e+(/.) - r K{\ + /.)r(/^) , 

c{\) = - r K{\ + /i)e+(/i) - r K{\- fi)rif-i) . (3.10) 

JAi JAi 



( P.SD is the generic form of a low energy spectrum in a system with two different branches of massless 



excitations. A finite size spectrum of the form ( [3.8D arises in many one dimensional systems, e.g 



integrable spin chains [ P2| , ^ [T^ and various models of correlated electrons [ P3| , P5| ] where the two 
critical degrees of freedom are holon and spinon excitations. In the present model it can be interpreted 
as a realization of two c = 1 Gaussian models. 

Decreasing the magnetic field further the universality class changes again at h = hd where Ai = 
and t>i vanishes leading to a divergence of the low temperature specific heat C ~ (ttT/G) (f f ^ + f 2"^) . 
At this point one of the massless modes has a quadratic dispersion and the system undergoes a 
Pokrovsky-Talapov transition into the U(l) Gaussian phase with central charge c = 1 that was al- 
ready identified for n < \ above. The value of h^ has to be determined by numerical solution of the 



integral equations (p7[). The complete magnetic phase diagram is shown in Fig. ^ 



This phase transition leads to a discontinuous change in the spectrum of critical exponents de- 
termining the long distance asymptotics of the correlation functions of the system. As an example 
we consider the spin spin correlator C^^{x) = {S^Sq). The massless magnon excitations discussed 
above lead to algebraically decaying correlation functions, from 03. 4| ) and ( P.9D we find that the lead- 



ing term in C^^ beyond the constants is of the form (1/x)" cos(2A;oa^) where /cq is a magnetic field 
dependent wave number and the exponent a is given by 

{2^2 for < /i< Ki 

E^mr + m') iorKi<h<K, 
For h -^ {hc2) one obtains a = 1 (2) from the appropriate expression as in region (i). In Fig. ^ 
the values of the exponent a as hd is approached from above and below are presented as a function 



of K. For i < K < 1 one observes a discontinuous increase of a at the transition h = hd. Such 
a change of the asymptotic behaviour of C^^ can be observed e.g. in the temperature dependence of 
NMR longitudinal relaxation rate, 1/Ti oc r""^ 

Finally we want to note that that these pronounced features in the /i-dependence of the magneti- 
zation and the critical behaviour of the system are reflected in the chirality properties of the system: 
Numerical computation of the chirality in the ground state as a function of the applied magnetic field 
for systems with up to 24 spins shows a smooth dependence of {x) on the magnetization. At the 
same time the fluctuations ((Ax) ) scale like 1/A^ for h < hd and are enhanced as the critical field 
is approached from below. Above hd they become smaller vanishing for the ferromagnetic state (see 
Figure ^. 

4 Summary and conclusion 

We have present results on the ground state properties and the magnetic phase diagram of a quasi 
one-dimensional spin-system with explicitly P and T breaking terms in the Hamiltonian. These 
terms generate chiral order in the system which in turn is found to lead to an interesting behaviour of 
the system when exposed to an external magnetic field. 

The question remains to which extent the observed behaviour is determined by the integrability of 



the system (|1.5D, (|1.6D. While a complete analyisis of this question is difficult we have diagonalized 
small systems to give a partial answer: 

Choosing the parameters as in ( |2.9[ ) with J2 = | the system interpolates between two soluble 



points: For k = the model reduces to the MG-model, k = ^ is the Bethe Ansatz soluble point. 
While the former one has a gap for spin excitations, the latter supports massless magnons. The van- 
ishing of this spin gap as a function of k has been studied using Lanzcos procedures and finite size 
interpolation. At the MG-point the gap above the two degenerate valence bond singlets is known very 
accurately, being Amg = 0.236 with finite size corrections scaling as A^ = A1^q{1 + const. /A^^) 
[P^. At K = I the gap vanishes as A ~ const /A^ according to (P^. It is difficult to find an interpo- 
lating expression between these limiting cases that leads to uniformly good fits of the numerical data. 
Still one can conclude that inclusion of the chiral term results in a rapidly vanishing gap, as can be 
seen in Figure ^. For k ^ 0.15 the numerical data do not allow to decide whether there is a finite gap 
or not. The ground state expectation value of the chirality shows a monotonic increase with k, with 
(x) ^ 0.58 at K = 1. 

Further studies are necessary for a better understanding of the intermediate phase transition in the 
integrable model at h = hd- ^ h = this transition occurs at the point k = 1, corresponding to a 



decoupling of the two sublattices. Whether such a 'dimensional reduction' is the origin of the phase 
transition at k < 1 and whether this feature can be used to describe the experimentally observed 
magnetic phases in the frustrated ABX3 compounds remains to be investigated in the context of 
systems with a larger number of coupled chains and eventually truly two-dimensional lattices. 

Acknowledgements 

We gratefully acknowledge useful discussions with H.-U. Everts, C. Lhuillier, U. Neugebauer, C. 
Waldtmann and A. A. Zvyagin. This work has been supported by the Deutsche Forschungsgemein- 
schaft under Grant No. Fr 737/2-2. The numerical calculations have been performed partly at the 
Regionales Rechenzentrum fiir Niedersachsen, Hannover, and the Zuse Rechenzentrum, Berlin. 

References 

[1] X. G. Wen, F Wilczek, and A. Zee, Phys. Rev. B 39, 11413 (1989). 

[2] M. Imada, J. Phys. Soc. Japan 58, 2650 (1989); H. Shiba and M. Ogata, J. Phys. Soc. Japan 
59, 2971 (1990); D. Poilblanc, E. Gagliano, S. Bacci, and E. Dagotto, Phys. Rev. B 43, 10970 
(1991). 

[3] S. Miyashita and H. Shiba, J. Phys. Soc. Japan 53, 1 145 (1984); S. Fujiki and D. D. Betts, Can. 
J. Phys. 65, 76 (1987). 

[4] T. Nikuni and H. Shiba, J. Phys. Soc. Japan 62, 3268 (1993). 

[5] S. V. Maleyev, Phys. Rev. Lett. 75, 4682 (1995). 

[6] Magnetic systems with competing interactions (Frustrated spin systems), edited by H. T. Diep 
(World Scientific, Singapore, 1994). 

[7] F D. M. Haldane and D. P Arovas, Phys. Rev. B 52, 4223 (1995). 

[8] C. Lhuillier, private communication. 

[9] M. Mino, K. Ubukata, T. Bokui, M. Aral, H. Tanaka, and M. Motokawa, Physica B 201, 213 
(1994). 

[10] N. StiiBer, U. Schotte, K. D. Schotte, and X. Hu, Physica B 213, 164 (1995). 



in 



[11] H. Bethe, Z. Phys. 71, 205 (1931). 

[12] C. K. Majumdar and D. K. Ghosh, J. Math. Phys. 10, 1388 (1969); ibid. 1399 (1969). 

[13] H. Godfrin, R. R. Ruel, and D. D. Osheroff, Phys. Rev. Lett. 60, 305 (1988). 

[14] H. Frahm, J. Phys. A 25, 1417 (1992). 

[15] V. Yu. Popkov and A. A. Zvyagin, Phys. Lett. A 175, 295 (1993). 

[16] C. Rodenbeck, Diploma thesis, Universitat Hannover, 1994. 

[17] H. Frahm and C. Rodenbeck, Europhys. Lett. 33, 47 (1996). 

[18] M. Fujii, S. Fujimoto, and N. Kawakami, J. Phys. Soc. Japan 65, 2381 (1996). 

[19] H. J. de Vega and F Woynarovich, J. Phys. A 25, 4499 (1992); B.-D. Dorfel and S. MeiBner, J. 
Phys. A 29, 6471(1996). 

[20] L. D. Faddeev and L. A. Takhtajan, J. Sov. Math. 24, 241 (1984), [Zap. Nauch. Semin. LOMI 
109,134(1981)]. 

[21] H. J. de Vega and F. Woynarowich, Nucl. Phys. B 251, 439 (1985). 

[22] A. G. Izergin, V. E. Korepin, and N. Yu. Reshetikhin, J. Phys. A 22, 2615 (1989). 

[23] F Woynarovich, J. Phys. A 22, 4243 (1989). 

[24] H. Frahm and N.-C. Yu, J. Phys. A 23, 2115 (1990). 

[25] H. Frahm and V. E. Korepin, Phys. Rev. B 42, 10533 (1990); ibid. B 43, 5653 (1991); N. 
Kawakami and S.-K. Yang, J. Phys. Condens. Matter 3, 5983 (1991). 

[26] T. Tonegawa and I. Harada, J. Phys. Soc. Japan 56, 2153 (1987); K. Sano and K. Takano, J. 
Phys. Soc. Japan 62, 2809 (1993). 



11 



Figure 1: Lattice on which the spin Hamiltonian ( |1 .3\ ) is defined: The two-spin exchange coupling is 
Ji and J2 on full and dashed lines, respectively. The three spin exchange oc xi,2 couples the spins on 
the corners of each triangle. 




Figure 2: Spinon dispersion of the integrable model (|1.6| ) for several values of k. 
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Figure 3: Chirality of the integrable Model (|1.6[). The system sizes are 2N = 8, 12, 16, 20, 24, cx3 
from top to bottom. 
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Figure 4: Magnetic phase diagram of the integrable system. Below hd the system has central charge 
c = 1, for hci < h < hc2 the model is a realization of two c = 1 models, and for h > hc2 the system 
saturates in a fully polarized (ferromagnetic) state. 
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Figure 5: Magnetisation of the integrable model vs. magnetic field for three values of k. 
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Figure 6: Filling fo the Fermi-Sea in presence of a magnetic field. Fig. (a): k < 1/4 and < /i < 
hc2- Fig. (b): k > 1/4 and h = hd, the position of the 'kink' in Figure |5|. 
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Figure 7: Critical exponent a of (S^Sq) as a function of k for magnetic field h 
(dashed line) and above (full line). 
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Figure 8: Chiral fluctuations (( A^) ) in the ground state for given magnetization a of the integrable 
system with k = | and n spins. The dashed line indicates the value of the magnetization at the critical 
magnetic field hd for the infinite system. 
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Figure 9: Spin gap A for Ji = 2 J2 = 1 as a function of the strength of the chiral field. The line is 
guide for the eyes and is given by a Gaussian. 
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